Sparse Exploratory Factor Analysis

Sparse principal component analysis is a very active research area in the last decade. It produces component loadings with many zero entries which facilitates their interpretation and helps avoid redundant variables. The classic factor analysis is another popular dimension reduction technique which shares similar interpretation problems and could greatly benefit from sparse solutions. Unfortunately, there are very few works considering sparse versions of the classic factor analysis. Our goal is to contribute further in this direction. We revisit the most popular procedures for exploratory factor analysis, maximum likelihood and least squares. Sparse factor loadings are obtained for them by, first, adopting a special reparameterization and, second, by introducing additional [Formula: see text]-norm penalties into the standard factor analysis problems. As a result, we propose sparse versions of the major factor analysis procedures. We illustrate the developed algorithms on well-known psychometric problems. Our sparse solutions are critically compared to ones obtained by other existing methods

Random forests highlight the combined effect of environmental heavy metals exposure and genetic damages for cardiovascular diseases

Heavy metals are a dangerous source of pollution due to their toxicity, permanence in the environment and chemical nature. It is well known that long-term exposure to heavy metals is related to several chronic degenerative diseases (cardiovascular diseases, neoplasms, neurodegenerative syndromes, etc.). In this work, we propose a machine learning framework to evaluate the severity of cardiovascular diseases (CVD) from Human scalp hair analysis (HSHA) tests and genetic analysis and identify a small group of these clinical features mostly associated with the CVD risk. Using a private dataset provided by the DD Clinic foundation in Caserta, Italy, we cross-validated the classification performance of a Random Forests model with 90 subjects affected by CVD. The proposed model reached an AUC of 0.78 ± 0.01 on a three class classification problem. The robustness of the predictions was assessed by comparison with different cross-validation schemes and two state-ofthe-art classifiers, such as Artificial Neural Network and General Linear Model. Thus, is the first work that studies, through a machine learning approach, the tight link between CVD severity, heavy metal concentrations and SNPs. Then, the selected features appear highly correlated with the CVD phenotype, and they could represent targets for future CVD therapies

Copper accumulation in vineyard soils: Rhizosphere processes and agronomic practices to limit its toxicity.

Viticulture represents an important agricultural practice in many countries worldwide. Yet, the continuous use of fungicides has caused copper (Cu) accumulation in soils, which represent a major environmental and toxicological concern. Despite being an important micronutrient, Cu can be a potential toxicant at high concentrations since it may cause morphological, anatomical and physiological changes in plants, decreasing both food productivity and quality. Rhizosphere processes can, however, actively control the uptake and translocation of Cu in plants. In particular, root exudates affecting the chemical, physical and biological characteristics of the rhizosphere, might reduce the availability of Cu in the soil and hence its absorption. In addition, this review will aim at discussing the advantages and disadvantages of agronomic practices, such as liming, the use of pesticides, the application of organic matter, biochar and coal fly ashes, the inoculation with bacteria and/or mycorrhizal fungi and the intercropping, in alleviating Cu toxicity symptoms

Piecewise smooth systems near a co-dimension 2 discontinuity manifold: can one say what should happen?

We consider a piecewise smooth system in the neighborhood of a co-dimension 2 discontinuity manifold $\Sigma$. Within the class of Filippov solutions, if $\Sigma$ is attractive, one should expect solution trajectories to slide on $\Sigma$. It is well known, however, that the classical Filippov convexification methodology is ambiguous on $\Sigma$. The situation is further complicated by the possibility that, regardless of how sliding on $\Sigma$ is taking place, during sliding motion a trajectory encounters so-called generic first order exit points, where $\Sigma$ ceases to be attractive. In this work, we attempt to understand what behavior one should expect of a solution trajectory near $\Sigma$ when $\Sigma$ is attractive, what to expect when $\Sigma$ ceases to be attractive (at least, at generic exit points), and finally we also contrast and compare the behavior of some regularizations proposed in the literature. Through analysis and experiments we will confirm some known facts, and provide some important insight: (i) when $\Sigma$ is attractive, a solution trajectory indeed does remain near $\Sigma$, viz. sliding on $\Sigma$ is an appropriate idealization (of course, in general, one cannot predict which sliding vector field should be selected); (ii) when $\Sigma$ loses attractivity (at first order exit conditions), a typical solution trajectory leaves a neighborhood of $\Sigma$; (iii) there is no obvious way to regularize the system so that the regularized trajectory will remain near $\Sigma$ as long as $\Sigma$ is attractive, and so that it will be leaving (a neighborhood of) $\Sigma$ when $\Sigma$ looses attractivity. We reach the above conclusions by considering exclusively the given piecewise smooth system, without superimposing any assumption on what kind of dynamics near $\Sigma$ (or sliding motion on $\Sigma$) should have been taking place.Comment: 19 figure

High Prevalence of Autoimmune Diabetes and Poor Glycaemic Control among Adults in Madagascar: A Brief Report from a Humanitarian Health Campaign in Ambanja

Madagascar is a geographically isolated country considered a biodiversity hotspot with unique genomics. Both the low-income and the geographical isolation represent risk factors for the development of diabetes. During a humanitarian health campaign conducted in Ambanja, a rural city in the northern part of Madagascar, we identified 42 adult subjects with diabetes and compared their features to 24 randomly enrolled healthy controls. 42.9% (n = 18) of diabetic subjects showed HbA1c values ≥ 9.0%. Unexpectedly, waist circumference and BMI were similar in people with diabetes and controls. Different from the healthy controls, diabetic subjects showed a low prevalence of obesity (5.7% versus 30%, p = 0.02). Accordingly, we found a high prevalence of autoimmune diabetes as 12% of people with diabetes showed positivity for the autoantibody against glutamic acid decarboxylase. Diabetic subjects with positive autoantibody had higher HbA1c values (11.3 ± 4.1% versus 8.3 ± 2.6%, p = 0.03) compared to diabetic subjects with negative autoantibody. In conclusion, here we describe the presence of diabetes and its features in a rural area of Northern Madagascar, documenting poor glycaemic control and a high prevalence of autoimmune diabetes. These data highlight that the diabetes epidemic involves every corner of the world possibly with different patterns and features

Early motion and directed exercise (EMADE) versus usual care post ankle fracture fixation: study protocol for a pragmatic randomised controlled trial

Background: Following surgical fixation of ankle fractures, the traditional management has included immobilisation for 6 weeks in a below-knee cast. However, this can lead to disuse atrophy of the affected leg and joint stiffness. While early rehabilitation from 2 weeks post surgery is viewed as safe, controversy remains regarding its benefits. We will compare the effectiveness of early motion and directed exercise (EMADE) ankle rehabilitation, against usual care, i.e. 6 weeks’ immobilisation in a below-knee cast. Method/design: We have designed a pragmatic randomised controlled trial (p-RCT) to compare the EMADE intervention against usual care. We will recruit 144 independently living adult participants, absent of tissue-healing comorbidities, who have undergone surgical stabilisation of isolated Weber B ankle fractures. The EMADE intervention consists of a non-weight-bearing progressive home exercise programme, complemented with manual therapy and education. Usual care consists of immobilisation in a non-weight-bearing below-knee cast. The intervention period is between week 2 and week 6 post surgery. The primary outcome is the Olerud and Molander Ankle Score (OMAS) patient-reported outcome measure (PROM) at 12 weeks post surgery. Secondary PROMs include the EQ-5D-5 L questionnaire, return to work and return to driving, with objective outcomes including ankle range of motion. Analysis will be on an intention-to-treat basis. An economic evaluation will be included. Discussion: The EMADE intervention is a package of care designed to address the detrimental effects of disuse atrophy and joint stiffness. An advantage of the OMAS is the potential of meta-analysis with other designs. Within the economic evaluation, the cost-utility analysis, may be used by commissioners, while the use of patient-relevant outcomes, such as return to work and driving, will ensure that the study remains pertinent to patients and their families. As it is being conducted in the clinical environment, this p-RCT has high external validity. Accordingly, if significant clinical benefits and cost-effectiveness are demonstrated, EMADE should become a worthwhile treatment option. A larger-scale, multicentre trial may be required to influence national guidelines. Trial registration: ISRCTN, ID: ISRCTN11212729. Registered retrospectively on 20 March 2017

A Study of Time-Dependent CP-Violating Asymmetries and Flavor Oscillations in Neutral B Decays at the Upsilon(4S)

We present a measurement of time-dependent CP-violating asymmetries in neutral B meson decays collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at the Stanford Linear Accelerator Center. The data sample consists of 29.7 ${\rm fb}^{-1}$ recorded at the $\Upsilon(4S)$ resonance and 3.9 ${\rm fb}^{-1}$ off-resonance. One of the neutral B mesons, which are produced in pairs at the $\Upsilon(4S)$, is fully reconstructed in the CP decay modes $J/\psi K^0_S$, $\psi(2S) K^0_S$, $\chi_{c1} K^0_S$, $J/\psi K^{*0}$ ($K^{*0}\to K^0_S\pi^0$) and $J/\psi K^0_L$, or in flavor-eigenstate modes involving $D^{(*)}\pi/\rho/a_1$ and $J/\psi K^{*0}$ ($K^{*0}\to K^+\pi^-$). The flavor of the other neutral B meson is tagged at the time of its decay, mainly with the charge of identified leptons and kaons. The proper time elapsed between the decays is determined by measuring the distance between the decay vertices. A maximum-likelihood fit to this flavor eigenstate sample finds $\Delta m_d = 0.516\pm 0.016 {\rm (stat)} \pm 0.010 {\rm (syst)} {\rm ps}^{-1}$. The value of the asymmetry amplitude $\sin2\beta$ is determined from a simultaneous maximum-likelihood fit to the time-difference distribution of the flavor-eigenstate sample and about 642 tagged $B^0$ decays in the CP-eigenstate modes. We find $\sin2\beta=0.59\pm 0.14 {\rm (stat)} \pm 0.05 {\rm (syst)}$, demonstrating that CP violation exists in the neutral B meson system. (abridged)Comment: 58 pages, 35 figures, submitted to Physical Review

Measurement of the Branching Fraction for B- --> D0 K*-

We present a measurement of the branching fraction for the decay B- --> D0 K*- using a sample of approximately 86 million BBbar pairs collected by the BaBar detector from e+e- collisions near the Y(4S) resonance. The D0 is detected through its decays to K- pi+, K- pi+ pi0 and K- pi+ pi- pi+, and the K*- through its decay to K0S pi-. We measure the branching fraction to be B.F.(B- --> D0 K*-)= (6.3 +/- 0.7(stat.) +/- 0.5(syst.)) x 10^{-4}.Comment: 7 pages, 1 postscript figure, submitted to Phys. Rev. D (Rapid Communications